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ABSTKACT 

Spia-angular momenturn transfer (or $pin-traxisfer for short) describes the angular momentum exchdoge between 
a spin-polarized current and a ferromagnetic conductor. When the conductor dimensions are reduced to around 
XOOxun Or below, the spin-angular momentum transfer effect becomes significant compared to the current-induced 
magnetic field. This paper describes some recent spin-transfer experimental findings in sub-lOOnm current- 
perpendicular spin- valve systems consisting of Co-Cu-Co nauopillars. The spin-transfer current is shown to 
cause a magnetic reversal of the thinner magtietic layer inside the nanopiUar. The reversal is experimentaUy 
shown to reach sub-nanosecond speed. The effect of spin-transfer is best understood in terms of its modifioation 
to the effective Laadau-Lifshis-Gilbert damping coefficient, either increasing or decreasing its value depending 
on the direction and magnitude of the spin-polarized current. For sufficiently laxge spin-current, the net damping 
coefficient may change sign, resulting in amplification of magnetic precession, leading to a magnetic reversal. 
At finite temperatures, the effect of spin-transfer is to either increase or decrease the thermal agitation of the 
nanomagnet. A quantitati\'e model is developed that adequately describes the finite temperatme experimental 
observations of the dynamic spin-traaisfer effect. 

Keywords: spin-transfer, magnetic nanodtructures, spin-vahre, finite temperature 

1. DNTTRODUCTTON 

Spin-angular-momentmn-transfer^ ^ (or spin- transfer for short) induced magnetization reversa) is a relative)y 
new phenomenon only observed in magnetic stnictureg smaller than 0,1 ^mi in size,^"® It reveals the interaction 
between magnetization and spin-polariaed transport current. The spin-transfier mechanism id further unambigu- 
ously observed in recent reports of spin-polarized current induced microwave emission in nano-pillar spin^valve 
structures.^^*^^ This paper give a brief discussion of the behavior of this representative spin-transfer device: 
namely a |Co|Cu|Co| current-perpendicular (CPP) spio-valve nanopillar. 

2. FABRICATION METHODS 

There ere commonly used fabrication methods for such junction devices. One is the so-caUed subtractive 
process, where the relevant multilayered metal film is formed first via vacuum deposition. The pillajr-shaped 
device structure is then formed by selectively etching away portions of the film. Alternatively, the pillar-shaped 
structure can also be formed via o» additive process. Junctions used in this study were fabricated using either 
a subtractive method^* or a variation of an additive process, namely a batch-fabricated stencil hard- masked 
substrate approach. The batch fabricated stencil substrate process is used in our lab for the fabrication of current- 
perpendicular spin-valve structures. Some details of this process has already been published elsewhere.*^^*^* lb 
briefly repeat, the process is summarized by Figure 1. 

The stencil Up typically is made of 20 xun or so thick platinum- An example of the top stencil structure is 
shown in Figure 2(a). A typical completed stendl structure woyld consist of ||10 Ti|100 ChiH Ti|50 SiOa |1 Ti|20 
Pt (I (numbers arc layer thicknesses in nm). 
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(a) Ion mill to open metal oiask (b) Wet eech to open 



(c) Deposit multilayer filxm 



(d) Form top level contact 
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Meta] 




FMI 
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Pt/Cu cnca»l 



Fi^re 1. Schematic of the btczicil-inask process, (a) E-bcam lithograpl^y And pattern transfer into the Pt stencil; (b) 
Wet etch 10 open up the insulator spacer and creating a controlled amount of undercut; (c) Deposition of magnotic stack 
foUon'Bd by metallic llUing that gives top electrode contact; (d) Optical lithography for the defi/iition of wiring. Figure: 
reproduced from ref.^^ 

CPP junctiond fbr this stttdy have been fabricated using planar magnetron sputter depodition, A typical 
junction stack consists of ||BottOTn|10 Cu|3 CojlO Cuj 12 Co|200 Cu|10 Pti|» A cross-aectional tranarni$$iOw 
electron microscopy view of the pillar shaped device fonwed through a. stencil structure is shown in Fig.2(b). 
This cross-section sample is cut from aii actual device using focused ion beam. 

After deposition of the layers fomung the junction device, optical lithography is used to form the necessary 
bottom and top electrodes 60 that the current passes through the stendl-defiued pillar. For reliable measurements 
of the junction resistance, the junction resistance needs to be well above the sheet resistance of the electrodes. 
A bottom electrode made of lOOnm copper typically gives a sheet resistance of about 0.08 to 0.1 0, allowing 
quantitative measurements of junction resistances down to about 1 



Quasi-static trajisport meaSure^nents were taken in a 4-probe geometry. A dc bias current is applied to the 
junction. At the same time a email ac current, ai)out lOO;iA rms at 331 Hz, is superimposed. The ac voltage 
is read out through a lock-in setup which measures the dynamic resistance of the junction at the given dc 
bias current. In what follows unless otherwise specified, all measured junction resistances refer to its dynamic 
resistance measured using this setup. In a typical quasi-static measurement, the bias current is stepped at a rate 
around 0 2 to 2 mA/min, while the magnetic field sweep rate (if swept) is on the order of 500-1000 Oe/min, 

A t3rpical ciurent-induced magnetic switching loop for a 0.05 x 0.15 ^tm^ junction is shown in figure 3(a). 
Genersdly speaking, the resistance is hysteretic both in applied field and in bias-current. The values of the 
switching threshold current and are a function of applied magnetic field and the sample's niagnetic 

and bias-current history. 

Figure 4 describes the hysteretic current and field dependence of the junction resistance. Here (a)-(d) are four 
directions when the bias-current is continuously swept in botli directions, while the magnetic field was stepped, 




Figure 2. (a)Scannbig electron microscopy view of the stencil mask, (b) A transmission electron microscopy cross-sectlona) 
image of a jHllar device deposited uinng planar magnetron i^uttering. 



3. QUASI-STATIC TRANSPORT BEHAVIORS 
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Figure 3, Left: An R {B) plot of a 
0.05 X 0.15 ^im^ jimctioR. Right: A 
A(/) pU>t. Inset shows the arran^ 
ment of the leads. Positive direo- 
tion of current is defined as current 
flowing from the bottom electrode 
thorough the pillar into the top elec- 
trode, dectrons fiow in the opposite 
direction. Measwement was done at 
ambient temperature. The vertical 
dashed line in the left-dlde panel in* 
dicates the bias: field at which data 
for R{I) yv&Tf: taken. 



either up or doiwii> between each bias current sweep. The direction of the current and field sweep is indicated 
by the arrows on the upper-right comer of each plot. Light color in the contour represerits hi^ resistance, dark 
color, low resistance. 

Tx^'o different regions can be seen from the data presented in Figure 4. In the low-magnetic field region, the 
/?(/) sweeps are hysteretic, with and Rildown) tracing out ft different boundary (subscripts refer to the 

direction of current sweep). In higher fields, the two brajiches of boundariea defined by R(I^) and R{f(tmtm) 
converge into a set of reversible bright or dark streaks. The^e strealcs correspond to dynamic resistance peaks or 
vialleys, caused by non-hysteretic steps up or down in the qua3i-statjc / - V diaracteristics. Notice the reversible 
Steps only appear in this case in the negative current direction, corresponding to electrons flowing from the thicker 
oobalt layer into the thimier. When one reverses the order of the thick and thin cobalt layer, the bias-current 
for these steps changes sign. 

The simple zero-temperature mono-domain model as described below for spin-transfer predicts a switching 
current threshold (albeit a threshold only for instability in the linearized solution of the Landau-Li&hit2-Gilbert 
equation) of 

Tc (H) = (aWMO {H^fff,^ 2nM,) (1) 

where the applied field H is along the easy-asds direction as defined by H^' a and b are the lateral Si^ses of Che 



f 




Figure 4. Contour plots of switching 
bouhdaries of a 0.05 x 0.10 fon^ junc- 
tion. Ambient temperature, field ap- 
plied along the eaiiy axis direction. The 
current was swept one full circle at a 
constant bi^ ficsld. The b'w field waa 
then stepped to the next vidue. The ar- 
rows at the uppuT right comer of each 
panel indicate the ciuTent sweeping and 
field stepping direction. The daslied 
line in (a) indicates the portion of data 
ivhere the slope-to^intercept ratio were 
taken for the field-dependent switching 
current /c(/f). 
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nanomagnetj t its thidoiess, M, its saturation magnetization. 7} is the spin-polarization of the current. Equ.(l) 
suggests the following relationship in the threshold current's magnetic field dependdtce: 



Namely* it predicts a linear dependence of on applied field, and the intercept to slope ratio (ISR) of the linear 
dependence should give Hk + 23rMa. For cobalt, Ms « 1, 400 emu/cm^ at ambient temperature,^^ which would 
predict an ISR of the order of 9 x 10* Oe even if <s: 25rifcf,. 

A nearly linear dependence of 7c as a Function of applied field H is indeed obeexved in a cetrtain range of Hy 
as shown in Figure 4 (a) -(d), an example is shown in Figure 4(a) by the dashed line. Tlus particular /dH) 
dependence gives an ISR of only 430 Oe, which i8 well below the value of close to 9»000 Oe aa expected from 
Eqn.(l)» A set of iSRs obtained from various samples is tabulated in Table 1. The experimentaJly observed 
values listed in Table 1 systematicaUy fall well below A^M^. In such measurements the uncertainty of zero- field 
H is no larger than the full demagneti^axion field from the piUar^s fixed cobalt layer, whose shape aoiiisottopy 
(which is the main cause of edge-coupling field) at 120 A by 0.03 x 0.15 fxvn? is no more than 2.5 kOe (Calculated 
using, for example, Eqn.(2) in Sun, ei oi.^^), which is insufficient to account for the discrepancy. It will be shown 
below that a large contributor to this discrepancy is the effect of fiuiite temperature. 

In addition to the simple linear dependence of /c On there are regiOJus in the (E, /) plane where non- 
hysteric steps were present. Such steps, when measured using the ac^diffeiential method give rise to the sharp 
peak (or valley, depending on the direction of the step) structure in the observed R{Hi /), appearing as bright 
or dark streaks in the contour plot. The detail of these peak width and height depends on the details of sample 
condition and the ac excitation current amplitude.'-^ A detailed analysis of the step width and related peak 
width and peak height wa$ recently given by Urashdin et ot^^ Alternatively, experimental evidences also exist 
that suggests the involvement of non-zero-k value spinewave excitations, with the wavelength of the excited 
spin-wave mode being simple multiples of the lateral sixe of the junction pillar.^* Details of magnetic excitation 
in this reversible region is likely to be complex, and T,vill depend sensitively on the exact nature of the magnetic 
boundary conditions of the pillar Junction structure. 



4. SWITCHING SPEED MEASUREMENT 

Direct measurement of the switching speed of spin-transfer junctions is no n- trivial because of the small signal 
level involved. For current-perpendicular (CPP) spin-v^xlv^ even with lateral sizes reduced to lOOmii or below, 
the junction resistance sill is only 1 — 50, and the raaguetoresistancc change is even smaller - usually only 
about about 3 ^ S % of total junction resistance. This results in a MR-related voltage signal typically on the 
order of 0.1 mV. Dynamic csJculations^^ give the generic time-scale of the reversal to be around (2?rA/j,) 7. where 
7 = 2pb/^ ^ gyro-raaen«ti<5 ratio, and M, is the saturation magnetization of the thia film nanomagnet. 
This estimate places the switching time ox* the order of 1 10 ns. 



Tb.ble 1. A list of mceeured slope, intercept and ISR for junctions with different size- (t)For comparison, one data^point 



Size (M^n^) 


Slope (AyOe) 


Intercept (mA) 


ISR (Oe) 


0.07 X 0.14 


6.46 x 10-^ 


7.34 


1,136 


0.07 X 0.14 


6.3^ X 10-^ 


5.87 


918.6 


0.05 X O.M 


7.81 X 10-^ 


3.83 


490.4 


0.05 X 0.10 


8.63 X 10"® 


7.14 


827,3 


0.05 X 0.10 


2.82 X 10-5 


7.81 


2,769.5 


0,05 X 0.10 {300K) 


1.31 X 10-^ 


5.63 


429-8 


0,05 X 0.10 (13-4K) 


1.08 X lO-*^ 


13.6 


1,255 


0.08 x 6.16 (360K) 


1.96 X 10-' 


17.6 


898.2 


0.08 X 0.16 (13.4K) 


5.9S X 10-^ 


10.8 


1,^9 


0.06 x 0.12(t) 


2.62 X IQ-** 


3.44 


1,314 
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Figure 5. (a) Quasj-statlc differentUl ceeistsnceof the 
jxiBCtion used for s^ritehing- speed sseasurement. (b) 
Quasi'Statie I-V character iHtics of the junction, with 
its linear (current-Independent) resistance subtracted 
out. Letters indicate the pulse sequence for higU-9peed 
switching measurement, m explained in the text. Fig- 
ure reproduced from Koch ei al.^^ 



.6 1 6 

Dvssf ou)n B* 

We focus on switchiAg-speed measurements performed at ambSent temperature. Pig.5(a) dhows the dynamic 
resSstancG of a typ!caJ sample used for switdiing-speed measurements. This particular junction presented here 
was patterned using a subtractive method.^ The active film stack sequence for this device is |Co|Cu|Co|. The 
free layer Co is about 3nm thick, and the lateral junction size is designed to be 0-05x0-10 fjJtn^. The average 
resistance for this junction was 3.7 Ohms and the MR resistaxice cliange is e-bout i percent of the total a small 
signal for broad-band measurement with severe common-modo j:oterfC5rence hoax the luain excitation current 
pulse. Additionally, samples often exhibit small changes of junction resistance over the course of a measurement 
session, perhaps due to contact condition changes. These add a few percent of long-term fluctuation to the 
measured resistance. To bring out a good broad-band MR signal in such a measurement environment requires a 
carefully designed multi-trace signal average process- 
Figure 5[b) plots the difference between the measured voltage and a straight line vs. the measured current. 
Detaib of the measurement has been published by Koch el^^ To briefly repeat, the measurement was done by 
alternating two pulse sequences, the •'signal'"' and the "background" sequence. The sample is dc biased at the 
current at points H (home) or B (background). For the "signal" sequence, a 1 microsecond pulse is ^plied to 
place the sample at point R (reset), fbr current-^dependence nieasiirements, a 300 (or sometimes a 30) ns pulse of 
varying amplitude is applied to place the sample at point P (pulse). If the pulse is sufficiently large, the sample 
will eventually switch to point S (switched). The voltage across the sample during this sequence is recorded using 
the TK6$00 oscilloscope. The background sequence is taken the same wy as the signal except without the reset 
pulse. When the sample is current pulsed nera the switching point, the lack of a reset pulse assures its moving 
from point 3 to point B during the backgrotmd sequence. The difference between the average voltages recorded 
between the two pulse sequences represents the magxvetimion reversal when the sample is biased between points 
P and S. Typically 10,000-traces are averaged in groups of 500 signal and background traces. 

Figure 6(a) plots the magnetization reversal time for a typical sample as a function of the pulsed current. 
'Pcom there a switching time r can be extracted. That is shown in Fig.6(b) and (c). Two poiirts are noted for 
the data shown in Pig.6(b-c)- First, at the high speed limit, the dependence of t"^ on bias current / is linear. 
Second, in the sub-threshold, largc-r regime, this linearity gives way to a curved onset which turns out to be 
exponentially dependent on the bias current as shown in Fig.6(c). We show below that the linear vs. / 
dependence stems from spin-transfer angular momentum conservation, and the curved onset relates to thermal 
activation. Both can be adequately described by the spin-tronsfer dyruunics with thermal noise. In addition, the 
threshold current observed in Fig.6(b) is about 1 1 mA for the threshold correeponding to the right-side (positive 
current) step of Fig. 5 which reads only about 4.5 mA. The difference between these two values lie in the vastly 
different time-scales of the two measurements- 
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Figure 6. (a)TimQ-dfip€ndent switching probability extracted from measurcmcntii illustrated in Fig.5; (b)Switchiiig speed 
t"* extracted from (a). Horizontal dashed lines \vith labeb indicate the switching ^eed corresponding to a rerversa] time 
of 1 and 2 nanosccondSj respectively; (c)Switching speed plotted on log-linear scale. Open and closed circle symbols 
represent the switching threshold on the positive and the negative current step of the junction shown in Fig.5. 

5, A MONO-DOMAIN MODEL FOR SPIN-TRANSFER MAGNETO-DYNAMICS 
5.x. Zero-tempetrature model 

Deitne a macro-spin with its magnetic moment m with a direcu)0« described by a \mit direction vector xim = 
»tn 01^) = sin e sin tpGj, -f- sin 6 C08 (pe,, -h cob ^e^ where S and a-te direction angles in a polar coordinate sy^enu 
The moment ra is situated in a combined energy potential oiV ^U{9,tp) that includes all energy-conserving 
torques m coiBtes to experience. The norraalijsed gradient of J7, expressed hi terms of ^cjj = (I/tti) VU 
includes terms applied magnetic field H, a uniaxial anisotropy whose strength can be characterised by a uniaxial 
aiusotropy field Hkt easy-plane anisotropy field that could be \xsed to describe a macrospin in thin^ 

film geometry, experiencing the demagnetization effect from the flat thin film geometry. Here the operator 
V = ep {d/dd) + (1/ sin 9) (B/dtp), with unit vectors eg and e^, denoting the direction of rot&tion for $ and ip, 
respectively. Note that other thj^ the applied field H, He// here 1$ in general not a dimple magnetic field vector 
but would be a fimction of the angular position of m. 

The procession dynamics of the macro-spin m uuder potential well U in the classical limit can be described 
by the phenomenolo^cal Landau-Lifehitz-Gilbert equation: 

= »«[H.„-0».«H.„] (3) 

where 7 2/ia/7t is the gyxo-maijnetic ratio. 

Assume a spin-polarized charge current / flowing along direction with spin-polarization along the imit- 
vector direction n, of the value v ^ {Ir - ^i) I (A + -^i) (where and 1^ designate current in spin channels 
parallel and auti-parallel to the direction of n,). This spin-polarized current is assumed to pass through m, 
beiiig re-polarized by m and hence acquiring a new 6pin-ix>iarization direction n^^ upon exiting. In this process 
a certain amount of spin-angular momentum from tlie spin-polarized current is absorbed by m. For realistically 
achievable current densities and magnetization' Ui metal-baaed magnetic canductor5> it is reasonable to assume 
the magnitude of m will remain unclianged during; the process, and hence to the lowest order, only the transverse 
component of spin-angular momentum deposited into m is ali^rbed by m, in the form of a precessiona] motion 
for xn. This is equivalent to adding a new term to £qn.(3], »o the precession dynamics of the macio-spln under 
the influence of a spin-transfer current will now read 

(i) ^ = - - (S) ™ X + «'>] 

where 5= h} (7i/2e] (l/ma) is the spin-angular momentimi transfer term. 

To understand the basic physics involved, examine the simple case when only an applied magnetic field H 
is present in He//, and that H and are collinear. That is, the incoming spin-current has its polarization 
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m X H 



(5) 



aligned (either parallel or ant]>parallel) to the direction of the applied field. In this simple limit one may write 

with 5 = Q + fti?// {2emHa) = a (1 H- I/Ic)y where 1^ — C^e/h) {A/t}) (abtMa) Hay a special case of the threshold 
current os deifined in Eqn.(J). This reveals the role of spin-polajized current /: it modifies the effectix'^e damping 
coefficient of the nanomagiJet. When the eflFective damping o becomes negative, the nanoma^et will amplify, 
rather than damp out, any disturbance away from itd equilibriuxn position, resulting in a magnetic instability 
leading to a magnetic re\'ersal. 

Fbr more general situations, since He// contains the an^ax position (B,ip) of m, a full stability analysis of 
Eqn.(4) is required. This in small cone-aagle limit < tt) can >e done ajnaJytically with a linearized Eqn.(4), 
as was done in ref.(^^). When averaged over a time-scale longer corapaied to the natural precession period 
Qk = 7^*,^^ it gives an effective damping coefBcient of 5 = q + hr}!/ [2em (J? + //t + 27rM5)| that describes 
the average cone- angle evolution {$ (t)). Here an uniaxial anisotropy field Hk and an orthogonal easy-plane 
enisotropj' term are included. A thin-film demagnetizatiou-reluted easy-plane anisotropy energy b assumed, 
therefore M9 = m/v where v is the volume of the nanomagneb. The resulting instability threshold, with a 



\miaxial anisotropy term and an orthogonal easy-plane anisotropjv' term included id 



Notice this threshold current Ic (H) depends linearly on applied field H. It should have a intercept-to-slopo ratio 
of Hfi 4- SttA/^. These can be compared with measurements. 

One set of experimentally measured ISR is listed in Table 1. It i 
ISK fieills well below the prediction of Eqn.(6). Furthermore, 
temperature of 14K appears to dhow less deviation than those : 
that finite temperature effect is at least partially responsible for 
result of Eqn.(6) and experimental observation. 



5.2. Role of finite temperature 

To model the finite temperature environment one follows the approach of Brown et al.^^ and Grinstein et al.^° 
by adding a Lflmgevin random field Hl to the efifective magnetic field term H^//. 
temperature T ss Hc^i - y/2akBT/ymIr^n,{ (0. = a;, y, z); where /, 



(6) 



is immediately apparent that the experimental 
(>ne observes that the ISR measured at a lew 
jasured at ambient temperature. This suggefita 
the discrepancy between the sero-tcmperature 



Hi relates to the system 
(£) is a Gaussian random function with 



the first two moments of {Iran (<)> = 0 and (i?an (t)) 1- The x^y^z components each has its own uncoxrelated 
Iran Without the spin-transfer effect the finite temperature LLG equation with a Langevin random field 
thus reads , ^ 

dm 



(7) 



This describes the dynamics of a macro-spin m sittinE in a potiutlal well U (fi,<p), with a thetmally 2ctiva,ted 
motion and a finite lifetime of staying inside the potential weJL namely a thermal life-thne r i^pioximately 
following the Boltzmann statistics of J 

where AI7 is the potential barrier height as seen from the local minimum where m fluctuates around, and 
TO ~ i/yHk 13 the reciprocal attempt frequency. At/ can be determined for any given potential landscape 
U(6,<p) and local minimum {^OjV^g). This approximation aliould work well when AU/koT '3> 1. The thermal 
activation life-time Eqn.(8) is detetxnined once a system is defin^ by Eqn.(7) and an energy landscape U (fl, v?)> 
represented in Eqn.(7) He// {9,<p) = (1/m) W 

The effect of finite temperature on the mOCro-Spin system's response to a spin-transfer excitation is two fold. 
First, it affects the average precession motion of the macro-spin 
a thermally distributed initial condition to the macro-spin. 



»y adding thermal activation. Secondly, it adds 
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5.3. ThermaJ activation under spin-transfer excitation 

Spin-transfer excitation adds an additiaaaJ torque, thus the zero temperature LLG Bon. (3) takes on the form of 
Ik3D(4). Similarly} after including the spin-transfer torque, and aesnme the spin-polarized current does oot cany 
significant entropy flow into the macro-Bpin system, one can re- write Eqn.(7) to read 

with H, representing the spin-transfer-torque related contribution^ as defined in £qn.(4). 

Solutions for Eqn(9) in general can only be obtained numerically. It ia however instructive to look at the 
special case when H, and H are both collinear to the easy-axis of the uniaxial anisotropy term Hj^ inside 
Hfi//, and the ea^y-plane anisotropy is zero. Id this case, and for ^mall cone-ang}e ^ <g: 1, one can show 
that He// = Hit + H, and hence the effect of is simply to modify He// — > H«// = {^jt ^ H •\- H9)ng = 
{Hi, + if) + 15^) = (l + T^tn) ^// ' re-write Eqn{9) as 

(;) ^ = - ■< i^"*^- - O (' * H7i«) » " H <"» 

Thus one may define a = (l -h jff^s) « = (1 + ^Ih) Qc, vnth, h = (2e/n) (oc/rj) {aJbtM^) (H + ff^), another 
special case of Eqn.(l),-and Eqn(XO) reverts to Eqn(7)*s form but with a renormalized by a: 

(^) ^ = " ^ + - (I) (") 

Other thajx the renormalSzation of o to 3, Eqn.(ll) is a mathematically equivaJent description of a macro- 
spin system as Eqn.(7) with the same amplitude of H^ — assuming the spin current is not by itself introducing 
additional thermal forces on the macro-spin^ which is reaj^onable as the amount of carriers involved in the 
paasing transport current per unit time is small comparoi to the equilibrium electron density. This then suggests 
f/L,i = Ht^i = y/ 2ak^T/ymIr(xn,i (*) remains valid for thi^ hypothetical macro-spui's LLG equation Eqxi(LO). 
However, in Eqn(lO) the damping coefficient o^is replaced by of. Tb maintain -ffL^i therefore it means the 
macro-spin would see an fictitious temperature T, such that af = oT. 

Since Eqn(lO) is equivalent to a macro-spin situated at temperature T with damping 5. one may further 
deduce that the thermal activation life time of the system can also be expressed in the form of Eqn.(8)» except 
with its temperature re-scalcd to T. That is: 

This conchision of a Uneaily current-dependent effective temperature caa be directly compared with experiment 
done in sub-tbreshoid driving current (|/| < 

5.4. Thermally distributed initial condition 

The initial condition of the macro-spin system under spin- transfer excitation is most important when the spin- 
transfer excitation is large - such as when it is near or exceeding the zero-temperature threshold current /c- in 
this case, when the spin-transfer current is applied suddenly at tone to, it will result in a fairly quick magnetic 
reversal. The switching time requited at 3!ero temperature is estimated to be:^^ 
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Here is the initial deviation of m from Jt$ es^^y-axis direction. The second line in Eqn-(13) is essentially a 
re^-^lation of angular monientum conservation. The first effect a iinite temperature brin^ is an initial deviation 
^0 ^th a thermal distribution. Thus the precise switching speed v/ill vary from meas\u:ement to measiircraent. 
At the same time, the thermal agitation during the course of reversal will add some uncertainty to the exact 
speed and trajectory of the reversal. This disturbance is liJkely to be small compared to the large oone-an{^le 
motion involved in these reversal events ati long as A^Z/it^T 3> 1 . 

To understand the consequerjice of a thermally distributed initial Oo one examines a special case with 

U {0, tp) = K (sin^ ^ + hp sin^ ^ cos2 - 2^ cos ^) (14) 

where K = m/fjb/2 is the uniaxial auisotropy energy constant, hp = AirM^Hk is the easy-plane anisotropy field 
in dimensionless unit, and h = H/Hk is the applied field, assuming a collinear geometry between H and Hfc. 
Assume further a small thermal fluctuation amplitude when i = K/IcbT » 1. In such limit, the main effect of 
finite temperature on delay r is through the initial axAgular position of 9: 



S^"" (?) ««P {-i^ (1 + h)) erf (TTtfVI^) dO 



(/-/«) 



(15) 



which in the limit of K/kgiT hp 1 gives a switchiiig speed based on the ensenible-avi^aged switching time 

(/-/c) (16) 



7Tln{K/kBT)]l /ip 

In this particular limit, switching speed has a current-dependence slope that is directly dictated by tbeawa) 
activation-induced initial angle. This ia only true if there's no transient disturbance fields during the application 
of the current pulse. Such transient fields would create a sudden rotation of the effective easy-axis dixectiou, and 
in effect creating a non-^ero and non-thermal initlaJ angle dictating the speed of spin-transfer switch. 

Summarizmg the temperature-dependence results discussed so far, we conclude: 

6. COMPARISON OF BXPBRIMBNTAL DATA WITH MONO-DOMAIN 

CALCULATIONS 

6.1. Spin-trdjfisfer^switching speed 

First, we examine the super-critical region. Experimental data from Fig. 6(b) gives for the positive side swltdiing 
a slope of dT~^fdI = 1.21 x 10^^ (sec • A)^^ According to Eqn.(17), assume spin polariaation factor 17 ^ 0.5, 
Mft = 1,400 emu/cm^ for cobalt, and a nano-magnet size of 0.05 x 0.10 /«n^ in the lateral direction and 30 A 
thick, havix)5 a magnetization m » 2. 1 x 10~^* emu, an uniaxial anisotropy field J^jb » 568. 6 Oe, and hence 
an anisotropy factor K = mHk/2 w 5.97 x 10"^^ erg which gives room temperature K/k^T a? 144, and an 
hp <l7rM«//f>. » 30.9. The exact value of applied field ^ in A » H/JIm is unknown in the experiment but is 
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of the order of Ht, and can be estimated to be /» 0.5. from the sub-threshold los-linear slope as is discussed 
belo..'. With these pamneters Eqn.(17) gives 6{rr' /dl ^ {^) ^^^^^j^^^^ « 1.3 iqu (sec-A)"!. 
This agrees with the experimental vahic shown in Fjg.6(b) to well within a factor of 2. 

Next we examine the sub-critical switching region. Details of this analysis is already published by Koch ei 
ai To briefly repeat, examine dat^ presented in Fig.6(c). The log-linear dependence of (r)"' on current is 
apparent. The slope of which according to Eqn.(17) gi^-es the strength of the pre^factor (1 - hf K/kuT, Mvhidh 
turns out to be 11. S7 and 28,78 for the uegative- and positive-side switching, respectively. 

Thus for the conditions tested experimentally, the functional farm of Eqn.(17) adequately describes both the 
i5ub- and the super-critical lixoit of the spin-transfer switching as a function of the switching current amplitude. 

6*2. Quasi-static switching boundaries and the role of finite temperature 

It can now be ahown that at finite temperatures, because of tliermal activation, the apparent measured quasi- 
static !c {H) will have its intercept- to-slope ratio fail below that of Hk-h2irM9, To see this recall the sub-critic*! 
expression from Bqn.(17). Witing A ^ K/kaT, the relaxation time at a given set of {HJ) hi the fT - / plane 
is 

r {H, I) - TO exp [a [1 - H/Hkf (1 J//c)] (18) 
Fbr an ensemble-averaged time-dependent magnetic moment m (i) this means 

'^0' «^ (1 - ff/f^k)^ (I - ///c)] (19) 

Fbr a time-dependent sweeping measurement of m (t), depending on the path of the sweep one can generally 
write 



^ (2) = ' exp [-A (1 - H/H.f (1 - ///J] (20) 

where v/v — dk/ dt and vj = dl/dt are (generally time-dependent) sweep rate for field and current respectively. 
The suffixes "0" and denotes the initial arid final point of the sweep. Eqn.(20) is a path integral with its 
value generally dependent on the path betweea (/io,ib) and (A/,//). Fbr example: for A » 1: 



^ { {^) acTtW [e-*<^-'^/^*>^f^-^/^»> ^ e-^<i-^/^0^] } , (/.^ep, constant H) 
<«P { (57k) [^e-^<^-^/^^)'(^-^/^.) - e-^<^-V/e)] } , (^.,^ep, constant 1) 



(21) 

For large A, the path-dependence is insignificant. The resulting contour plot of vrif/mo as a function of 
(ff,/) is illustrated by the numerical vesultfi in Figure T. From Fig. 7 it is apparent that at finite K/ksT ratio, 
the apparent switching boundary (as indicated by the dark to bright transition) bas a reduced slope-to-intercept 
ratio as those of zero-temperature Ic (/f), indicated in Pig.7(a) by the white line above the contour edge. This 
exercise demonstrates clearly the eflfect of finite temperature on the reduction of the apparent slope-to-hitercept 
ratio of the quasi-staticaDy measured Ic (H). This mechanism is at least partially responsible for the discrepancy 
between the zero-tempcrautre model prediction and the experimental observation shown in Table I. 

7. CONCLUSIONS 

Spin-traasfer-induced magnetic excitation and magnetic switching is manifested in sub-lOOnm current-perpendicular 
$pin-valves composed of |Co|Cu|Co| bilayers. The Spin-transfer excitation is slKwm to cause both hysteretic mag- 
netic switching aad non-hysteretic magnetic excitation. For magnetic reversal, finib-nanosecond switching speed 
is experimentally demonstrated. A finite temperature single-domain model is developed to successfully describe 
the observed magnetic dynamics. 
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Figure 7. Finite temperature switching probabil- 
ity contour as a function of current and field far 
two difiiorent valuer of K/k^T. Hdxe a cKmstant 
curreait Sweep rate of dl /dt 60 fiAJtsec is used. 
Thin contour plot does! not look significantly differ- 
ent for a field b-weeping approach with dH/dt ^ 20 
Oo/soc. 



I 0 



I 01 



8. ACKNOWLEDGMENT 



We wish to thank the valuable input from the IBM MRAM team. One of us (JZS) acknowledges a productive 
collaboration on topics related to this study with B. Oezylimaz and A. D. Kent at New York University supported 
by NSP-FRG Grant# 0101439- 



1. J. C. Slonc«€wski J, Magn. and Magn. Mat. 159, p. Ll, 1996. 

2. L. Berger Phys, Rev, B 54, p. 9353, 1996, 

3. J. C. Slonczewski J. Magn. Magn, MaUr, 195, p. L261, 1999. . 

4- M. Tsoi, A. G- M- Jansen, J. Bass. W.-C. Chiang. M. Seek, V. Tsoi, and P. Wyder Phys, Rev. Utt 80. 
p. 4281, 1998. 

5. J. Z. Sun J. Magn. Magn. Maicr, 202, p. 157, 1999. 

6. E. B. Myers, D. C. Ralph, J. A. Katine, R. N. Louie, and R. A. Buhrman Sdcncc 285, p. 867, 1999. 

7. J.-B. Wegrowe, D. Kelly, Y. Jaucatd, P. Guittienne, and J.-P. Ansermet Europhysics Letters 45, p. 626, 
1999. 

J. A. Katine, F. J. Albert, R. A. Buhrman, E. B. Myers, and D, C. Ralph Pkys- Rev. Lett- 84, p. 3149, 
2000. 

9. J. Grollier, V. Cros, A. Hawizic, J. M. George, H. Jaffr^, and A. Pert Appl Phys. Lett 78, p. 3663, 2001. 

10. S. I Kiselev, J. C. Sankey^ L. N. Krivorotov, N. C, Emley, R. J, Schoelkopf, R- A. Buhrman, and D. C, 
Ralph NcUure 425, p. 380, 2003. 

11. W. H. Rippftrd, M. R. Pufell, S. Kalca, S. E. Russek, and T. J. Silva arXiv:cond'mat/030904S , 2003. 

12. J- Z. Sun, J. C. Slonczewski, P. L. Trouilloud, D. Abraham, 1. Bacdius, W. J. Gallagher, J. Hummel, Y. Lu, 
G. Wright, S. S- p. Parkin, and R. H. Koch Appl Phys, Lett, 78, p. 4004, 2001. 

13. J. Z. Sun, D. J. Monraia, M. J, Rooks, and R. H. Koch Appi Phys. Lett. 81, p. 2202, 2002. 

14. J. Z. Sun, D. J. Monsrma, T. S. Kuan, M. J. Rooks, D. W. AbraJiacoi, B- Qezyilaz, A. D. Kent, and R. H. 
Koch J, Appl. Phys. 93, p. 6859, 2003. 

15. C. Kittel, p. 429, Introduction to Solid State Physics, 6th Ed. John Wiley k Sons, 1986. 

16. S. Urazhdin, N. O. Birge, W. P. P. Jr., and J. Ba68 arXiv:COnd-mat/0SlSS87v2 , 2003, 

17. J. Z. Sun Phys. Rev, B 62, p. 570, 2000. 

18. R. H. Koch, J- A. Katine, , and J. Z. Sun, to appear in Phy$, Rev. Lett. 2004. 

19. W. F. Brown Phys. Rev, 130, p. 1677, 1963. 

20. G. GrtTji$tein and E. H. Koch Phys, Rev, Lett, 90, p. 207201, 2003. 



PAGE 3(»30'RCVDAT9/11<20(I(11:5S:47PM [Eastern Daylight Tim^^ 



REFERENCES 



11 



